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bJO Abstract 

We present the first examples of black holes with only one Killing field. The solutions 
describe five dimensional AdS black holes with scalar hair. The black holes are neither 
stationary nor axisymmetric, but are invariant under a single Killing field which is tangent 
to the null generators of the horizon. Some of these solutions can be viewed as putting 
4^ black holes into rotating boson stars. Others are related to the endpoint of a superradiant 

"*T* instability. For given mass and angular momentum (within a certain range) several black 

hole solutions exist. 
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1 Introduction 

All known black hole solutions have at least two Killing fields. This is not just due to the 
fact that finding analytic solutions of Einstein's equation usually requires the assumption of 
symmetries, but also follows from theorems which say that stationary black holes must be 
axisymmetric [T1[21I3]- However, it has been proposed that black holes in anti-de Sitter (AdS) 
space might exist with just a single Killing field The argument is based on superradiance. 
Recall that a mode e""^*"'"*™''^ of a scalar field can increase its amplitude by scattering off a 
rotating black hole with angular velocity Qh satisfying u < rnQn- In asymptotically AdS 
spacetimes, this leads to an instability since the outgoing wave is reflected back onto the 
black hole and scatters again further increasing its amplitude. This process decreases Qh 
and eventually results in a black hole with "lumpy" scalar hair rotating around it. Such a black 
hole would not be stationary or axisymmetric, but invariant under just a single Killing field 
which co-rotates with the scalar hair. 
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The main goal of this paper is to numerically construct such black holes and study their 
propertied It will be convenient to work in five dimensions. This is because a clever ansatz was 
recently found for five dimensional scalar fields in which the coupled Einstein scalar field equa- 
tions reduce to ordinary differential equations even though the scalar fields are invariant only 
under one Killing field [H] • In [B] a potential was added for the scalar fields and asymptotically 
flat solutions without black holes were found. These were interpreted as rotating boson stars. 
We will set the potential to zero (so the scalars are massless) and add a negative cosmological 
constant. We will see that rotating boson stars still exist in this case. The first example of a 
rotating boson star in AdS was given in [7], where the authors constructed a cohomogeneity-1 
solution by working in three dimensions. 

Of more interest, we show that one can add small spherical black holes in the center of these 
boson stars without destroying the solution. There are results in the literature stating that one 
cannot add black holes inside boson stars [8l |9]. These results apply to static black holes and 
follow essentially because the boson star has e~"^* time dependence and t — )• oo at the horizon 
of a black hole. The scalar field thus oscillates infinitely often near the horizon and cannot be 
smoothly continued inside. We will see that if you put a rotating black hole inside a rotating 
boson star, this obstruction goes away. One way to see this is to note that a mode e~*'^*+*'"'^ is 
invariant under K = d/dt + {u / m)d / dip . If the black hole has angular velocity VLh = uj/m, the 
Killing field K is the null generator of the horizon, and the scalar field is essentially stationary 
in a co-rotating frame near the horizon. 

The special ansatz for the scalar fields corresponds to an m = 1 mode. For this mode we 
determine which Myers-Perry black holes are unstable to superradiance and find that they all 
have radii much less than the AdS radius. For the same energy E and angular momentum J, we 
find black holes with scalar hair which represent the endpoint of the superradiant instability for 
this mode. In the full theory with no restriction on the scalar fields, these hairy black holes are 
still subject to superradiant instabilities associated with higher m modes. In the final section 
we discuss the endpoint of these black holes in the full theory. 

Similar to the asymptotically fiat case, rotating boson stars in AdS have quite an intricate 
structure. As one increases the central density, E and J execute damped oscillations around 
critical values E^. and Jc- Thus, not only is there a maximum E and J, but there are many 
different boson stars with energy and angular momentum near Ec and Jc- For each of these 
boson stars, one can add black holes of various sizes. The area of the black hole can range 
from zero up to a maximum. This leads to a high level of black hole nonuniqueness. There can 
be several different hairy black holes with the same energy and angular momentum. We will 
see that as one increases the size of the hairy black hole, the solution approaches a singular, 
extremal configuration. 

2 Model 

2.1 Equations of motion 

We consider five-dimensional Einstein gravity with a negative cosmological constant minimally 
coupled to two complex massless scalar fields 11 , i.e. with action (we set Newton's constant 



^Another way to construct asymptotically AdS black holes with only one Killing field is to modify the 
boundary conditions at infinity so that the boundary metric is not conformal to a round sphere cross time [5]. 
In this case, the spatial symmetries are broken by the boundary conditions. We are interested only in black 
holes with the standard asymptotic behavior. 



3 



G = l) 



S 



IGtt 









2" 






vn 




Im 









(2.1) 



We will look for boson star and black hole solutions of this theory whose gravitational and 
scalar fields obey the ansatz: 



h[di, + ^d0 - ndt\ + ^ (d^2 + sin^ ^d02) 



n = Ue-''^'+'^ 
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, (2.2) 
(2.3) 



Here, f,g,h,Q and 11 are real functions only of the radial coordinate r, i.e. the ansatz is 
cohomogeneity-1. Surfaces of constant time t and radial r coordinates have the geometry of a 
homogeneously squashed S^, written as an 5*^ fibred over S"^ = CP^ . The fibre is parameterized 
by the coordinate ■?/', which has period 27r, while 9, are the standard coordinates on S'^. Since 
(f) also has period 27r, it might appear that 11 is not single valued, but the Hopf fibration requires 
that t/; —)■?/' + TT when — )■ + 27r. Thus 11 is indeed single valued on spacetime. 



The motivation for the particular ansatz (2.2) for gab will become clear later in this section. 
This metric has five linearly independent Killing vector fields, namely dt, and the three 
rotations of S"^. However, it is easy to check that the only linear combination which leaves H 
invariant is 

K = dt + ujd^. (2.4) 

Since a symmetry of the solution must leave both the metric and matter fields invariant, any 
solution with nonzero scalar fields will have only one Killing field. Of course, one does not 
usually expect there to be solutions in which the matter fields have much less symmetry than 
the metric. However the ansatz (2.3) is speciaj^in that its stress tensor 



Tab = dJi*dbTi + daTidbTl* - ^?ab 9,H9^H* 



(2.5) 



has the same symmetry as the metric (2.2). One way to see this is to think of the two complex 
components of H as coordinates on . The fact that |Hp is only a function of radius means 
that at each (r, t), H maps the squashed 3-sphere in spacetime into a round 3-sphere. The first 
term in the stress tensor is just the pull-back of the highly symmetric metric on this sphere 
and the second term is proportional to the spacetime metric. 

We thus expect our ansatz to have nontrivial solutions, and indeed the field equations 
Gah — G^~'^gab = Tab and V^H = 0, lead to a consistent set of ODE's: 
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(2.6b) 



^It was first considered in 
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^ i!:!^|p^ + 4^/,S = , (2.6f ) 

_ U\Q-uj)' h'ifh)' {r'hyz Sjfhr'Y ^ 2 -U' r'hQ'' _ 

r-ippK^g pg + 4/i2/ + r5/i/ 2/ir3/ ^ /£2 + r2/i/ 4/^ 

(2.6g) 

where = h + lP — jP' — 2, and ' denotes differentiation with respect to r. The last two 
equations can be regarded as constraint equations, which under the ffow of the remaining five, 
obey: 

C-f^ + ^l.a^.5V., (2.7a) 



j p h r 
■Eir^h)' /i2-n2 rh!{]h)' 3(fhr^y 2rh ,2 



Ci + ^fiVc2. (2.7b) 



h \r 



If we had considered a single complex scalar field, we could still assume a Fourier decom- 
position along the time t and azimuthal coordinates V^, of 5"^. This guarantees that Tab is 
independent of t, ip, 0. However the stress tensor would still have contributions depending on 
the polar coordinate 9. When plugged into the Einstein equation, such a Tab necessarily sources 
a gravitational field that is cohomogeneity-2, i.e. that has a radial and polar dependence. The 



simple ansatz (2.2) would then not be appropriate and finding a hairy solution would require 
solving a coupled system of PDEs; a much more difficult task. Note also that a more general 
ansatz for the doublet scalar field would consider a more general Fourier decomposition e™'^, 
i.e. with an arbitrary integer azimuthal quantum number m not necessarily restricted to the 
case m = 1 that we consider. Although we certainly do not rule out the possibility of finding 
a cohomogeneity-1 ansatz for the gravitational field also when m > 1, this is not a trivial task 
and we will not attempt it here. Moreover, as we will prove latter, the superradiant instability 
of Myers-Perry— AdS black holes is stronger for the m = 1 mode than for higher m modes. So 
the hairy (black hole) solutions with m = 1 that we will construct are the first to form from a 
generic perturbation. 

2.2 Boundary conditions for boson stars and black holes 



To discuss the boundary conditions of the solutions, we first notice that the ansatz (2.2) and 



(2.3) still has a residual gauge freedoom, that leaves the gravitational field (2.2) invariant: 



ip^ip + Xt, n^Q + X, u^uj-X, (2.8) 

for an arbitrary constant A. Unless otherwise statecQ we use this gauge freedom to require 
that f2 — )• as r — )■ 00. We are thus in a frame where the solution does not rotate at infinity. 



^In section |4] we will construct numerically the hairy solutions. We will find that in certain regions of the 
parameter space the numerical convergence is better if we work in the unphysical gauge where we choose A to 



5 



i.e. the conformal boundary of the hairy solutions is an Einstein universe R x S^. This also 
guarantees that, in the black hole case, the angular velocity measured at the horizon, Qh, is 
the appropriate quantity for a thermodynamic analysis [lOl [TTl [12] . 

The asymptotic boundary conditions are the same for the boson star and black hole. More 
concretely, we require that our solutions f,g,h,fl asymptote to global AdS with the next-to- 
leading order terms in the expansion (determined by constants Cf, Ch, Cn) being responsible 
for the generation of mass and angular momentum. The boundary condition for 11 is determined 
by the requirement of normalizability. This demands that the massless scalar field decays as 
r~^. In the above non-rotating frame, the asymptotic boundary conditions are then: 



»,2 nr p"^ r, 



h 



H™ = i + -^^ + ^(--')' ^Uoo = ^ + ^(--') . n|_ = — + o(r-^). 

(2.9) 

Note that e is a quantity that will be frequently used henceforth and is a dimensionless asymp- 
totic measure of the amplitude of the scalar field (condensate). 

The second set of boundary conditions naturally depends on whether we are looking for 
boson star or black hole solutions. Regular boson stars are smooth horizonless solutions (with 
harmonic time dependence), i.e. its boundary conditions at the origin are: 

/L^o = l+^(-')' ^L^o = ^(l)' H.^o = 1+^ ' ^L^o = ^(l)' n|^_^„ = 0(r). 

(2.10) 

Note that regularity of 11 requires 11 (r) to vanish linearly at the origin. 

For the black hole, the inner boundary is at its horizon where the function /(r) must vanish. 
We take this condition as our definition for the location of the black hole horizon r_|. (the largest 
root of /). The other functions must be regular at this hypersurface: 

I ^ ^~h) ' ^ I 7" ^T-\- ^ ) ' ^ I T° ^Tj^ ^"^^ ' I ^T°_[- ^"^^ ' I 7" J 

(2.11) 

Regularity of 11 on the horizon imposes a further restriction. Multiplying (2.6e) by and 



evaluating on the horizon shows that the frequency to of the scalar field must equal the angular 
velocity of the black hole Qh, i-^- 

u = Qh = (black hole case) . (2.12) 

Recall that the only Killing field of our solution is K = dt + ood^ (2.4). In general, this has 



norm \K\'^ = —fg + r'^h{uj — Q)"^. On the horizon, (2.11) and (2.12) imply K is null, so the 
event horizon is also a Killing horizon associated with the only Killing field of the solution. K 
is always timelike just outside the event horizon and in a neighboorhood of the origin r = for 



the boson star. The asymptotic boundary conditions (2.9) imply that li^l — )■ (cu^ — V^^) 
r — )■ oo. Therefore the Killing field of the hairy solution will be asymptotically timelike, null 
or spacelike depending on whether ui < 1 , ui = 1 or ui > 1, respectively. The solutions we 
find will all have ui > 1. Thus, the solutions are not globally stationary. There is an effective 
ergoregion at large radius. It turns out that the scalar field is actually concentrated in this 
region. 



be such that w = (which uTiphes il(oo) ^ 0). When we find it necessary to do the computations in this gauge, 
we will however present our final results in the physical gauge where U{oo) = (and thus a; 7^ 0). Note also 
that the boundary conditions in the gauge w = can be easily read from the boundary conditions for the gauge 



ri(oo) = once (2.8) is applied. 
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2.3 Asymptotic conserved quantities and thermodynamics 



Even though the vectors dt and are not Kilhng vectors of the full hairy solution, they are 
asymptotic Killing vectors (since the scalar field vanishes asymptotically). We can therefore 
define conserved quantities associated with these asymptotic symmetries. To compute the 
conserved quantities, we use the Astekhar-Das formalism [1^, which we now briefly review. 
Start by introducing a conformal metric ^ related to the physical metric g through 'cjab = 



A'^Qab- In our case, the physical metric is given by (2.2) and the appropriate conformal factor 
is A = 1/r. The conformal boundary surface Sa=o is defined by A = and has normal vector 
= daA. Let Cabcd bc the Weyl tensor associated to the conformal metric. Define the tensors 
Kahrd = A.^~'^Cabcd\jy_Q and Eab = ^'^Kabcd'n^n'^- On the conformal boundary Sa=o consider now a 



timelike surface with t =constant and with normal ta = dj:. This defines a (d — 2)-dimensional 
hypersurface with induced metric h that we call E. Take now an asymptotic conformal Killing 
vector S,. The conserved charge Qc associated to this symmetry is [13] 



d-3 



57r 



SabCt'dE 



(2.13) 



where ± for a timelike/spacelike conformal Killing vector, respectively. 

For solutions described by (2.2) and (2.3), we can take ^ = {dt, d^) and d = The 
associated conserved quantities are the energy, E = Qs^, and the angular momentum, J = Qa^, 
given by 



E 



IT . 



(4C, - 3C 



f) ' 



J 



(2.14) 



where the constants C/, Ch, Cq can be read from the asymptotic decay of the gravitational 



fields in the boundary conditions (2.9) 



Boson stars have no horizon and therefore they are zero entropy objects. They do not 
have an intrinsically defined temperature. If the solution is a black hole, we have just seen 
that the Killing vector K = dt + ud^ becomes null at the event horizon r = r+. This is 
therefore a Killing horizon and its temperature is given by the surface gravity divided by 27r, 



-iVKy\^ / The black hole entropy is, as usual, the area of the event horizon 

divided by 4. This gives 



H 



4:71 



S 



r+ 



(2.15) 



where ^4^.2 is the area of a unit {d — 2)-sphere. 

Any family of solutions to Einstein's equation with a Killing field satisfies an analog of the 
first law of thermodynamics. This follows from a Hamiltonian derivation of the first law HA 



When the Killing is given by (2.4), this takes the form: 



Boson star: dE = udJ , (2-16) 
Black hole: dE = udJ + TndS , with a; = fi// . (2.17) 



2.4 Myers-Perry— AdS solutions 



A very special black hole solution of (2.6) is the one where the scalar field vanishes. It describes 



a rotating black hole where the angular momenta along the two possible rotation planes in c? = 5 



7 



is the same. This solution is the d = 5 case of "Myers- Perry— AdS" (MP-AdS) black holes jlOj . 
with equal angular momenta and with cohomogeneity-1 that exist in odd dimensions d = 2A^+3 
[T5| [T6] . Although we are mainly interested in the d = 5 solutions, we find it useful to discuss 
briefly the higher- dimensional counterpart of these black holes. Therefore in this subsection, 
we present the results for an arbitrary odd dimension. The equal angular momenta MP-AdS 
black holes have a gravitational field given by 



(2.18) 



where 



/ = 1 + 



„2N 
M 
.2N 



a 



+ 



2Af 2 
' M " 

^2(Af+l) 



h = l + 



2Af 2 
^2(Af+l) 
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1 



n 



^2N 



(2.19) 

2{N+l)gab,SindA = AJx" 



and gab is the Fubini-Study metric on CP^ with Ricci tensor Rab 
is related to the Kahler form J by dA = 2J. 

Asymptotically, the solution approaches AdS space. The event horizon is located at r = r+ 
(the largest real root of /) and has a null tangent vector dt + ^nd^ where the angular velocity 
is: 

1 



n 



^2N 



H 



2N+2 , 



2iV^2 — ^ 



where 



Qext 



1 + 



;i + Ny^ 



(2.20) 



The solution saturating the bound in the angular velocity corresponds to an extreme black hole 
with a regular, but degenerate, horizon. Note that this upper bound in Qh^ is always greater 
than one and tends to the unit value in the limit of large r+ / i. 

It is convenient to parameterize the solution in terms of (r+, Qh) instead of {tm, a) through 
the relations, 

AN+l) 1^2 , .2^ ^2 ^2^ 

^ (2.21) 



^2N 

M 



a2 {rl + &) ' 



The temperature and the entropy as defined in (2.15) are given by 



H 



rl + £2 Nf _ (AT + l)rl {n^ f - 1) ^ ^ A^^^^TV+i 



2'nr^P 



P + ri{i-pn\) ' 

(2.22) 

and the temperature vanishes for the extreme configuration. The mass E and angular momen- 



tum J, defined in (2.13), are [12 



E 



A 



2N+1 2N 
871 "^'^ 



N + 



1 rxeni 



2{rl + &)\ 



J 



A 



2N+1 

Svr 



(2.23) 

In this 
0. The 

line element (|2.18|) for the equal angular momenta MP-AdS solution and the fact that it is 



The = 1 case describes the d = 5 black hole that we are most interested in. 
case ([2l8|)-([2l9| reduce to (|2g-(|23|) with CP^ = S^, 2A„dx" = cos^d0, and n(r) = 



superradiant unstable (thus indicating a possible bifurcation to a new branch of hairy black 



holes) motivated our choice of the gravitational ansatz (2.2) for the hairy solutions. Henceforth, 
unless otherwise stated, when we refer to the equal angular momenta Myers-Perry— AdS black 
hole we will assume that d = 5. 
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3 Perturbative results 



In this section we construct rotating boson stars and hairy black holes that have a single 
Killing field (2.4). This is a perturbative construction, i.e. valid only for small energy and 
angular momentum. 



3.1 Rotating boson stars 



Boson stars are smooth horizonless geometries with harmonic time dependence. There is a 
one-parameter family of solutions, which in the perturbative regime can be parametrized by e. 
Recall that this is a measure of the amplitude of the doublet scalar field at infinity (2.9). We 



can construct perturbatively the fields of these rotating stars through a power expansion in the 
condensate amplitude e around global AdS, 



F(r,e) = X^F, 



2i 



2i 



n(r, 



En: 



2j+l 



r e 



2j+l 



U2i e 



2i 



(3.1) 



where we use F to represent schematically one of the gravitational fields in (2.2). That is, each 
gravitational field F = {/, g, h, Q} has a power expansion in even powers of the condensate 
e, with expansion coefficients F2j{r) = {f2j{r), g2j{r), h2j{r),fl2j{r)} at order e^-'. The doublet 
scalar field has an expansion in odd powers of e. Note that we expand not only the gravitational 
and scalar fields but also the frequency. The reason being that at linear order the frequency is 
constrained to be one of the normal modes in AdS but, when non-linear effects are included, 
this frequency receives corrections. When e = 0, (3.1) describes global AdS spacetime. Using 



a standard perturbation theory strategy, we plug the expansions (3.1) in the equations of 
motion (2.6) to solve for the coefficients f2j, g2j, h2j,fl2j,^2j+i and U2j. At any order in the 



expansion, the fields must decay asymptotically according to the boundary conditions (2.9) 



and the solutions must be regular at the origin as dictated by the boundary conditions (2.10). 

The leading order contribution in the expansion, n = 0, describes the linear perturbation 
problem where we introduce a non-trivial massless doublet scalar in the AdS background, but 
this condensate does not back-react on the gravitational field. The most general massless 



doublet scalar field in AdS that decays with the power law given in (2.9) and satisfies (2.6e) is 
given by a sum of normal modes that are labelled by k and given by (2-F1 is the hypergeometric 
function). 



U(r) 



5/2 



5-ui 5 + uj£ 



uji = 5 + 2k, (fc = 0,l,2. 



(3.2) 

where the quantization condition on the allowed frequencies u follows from requiring the regu- 



larity condition (2.10) at the origin. This selects the discrete frequency spectrum of stationary 
scalar normal modes that can fit into the AdS box. The non-negative integer k describes the 
several possible radial modes or overtones. 



A non-linear boson star with Killing vector (2.4) is sourced by just one of the modes de- 



scribed in (3.2). In this paper, we will focus on just the lowest {k = 0) normal mode, since 
this is the one that describes the ground state boson star with the lowest mass to angular 
momentum ratio. For k = 0, (3.2) simplifies. In the nomenclature of (3.1), the full solution 



that satisfies the equations of motion (2.6) up to order 0{e) is then 



If we keep moving in the expansion ladder to the non-hnear level, the scalar field back- 
reacts in the metric (and vice- versa), with the lower order contributions sourcing the higher 



order fields. At each order n we have to solve the coupled system of equations (2.6) up to 



order O (e 



This can be done analytically. Typically, these are second order ODE's so the 
most general solution is a linear combination of two solutions for each field. The two sets of 



boundary conditions (2.9) and (2.10) fix the two coefficients of the linear combinations that 
yield the physical solution. At order and actually at any even order e^, the equation for 



the scalar field (2.6e) is trivially satisfied and the gravitational field equations (2.6a)-(2.6d) are 



the ones whose solution encodes the non-trivial information on the back-reaction of the scalar 
field on the gravitational field. On the other hand, the equations at any odd order e" are used 
to find the corrections to the scalar field. This construction can be carried with relative ease 
up to arbitrary higher orders in e with the aid of Mathematica. The expressions for the fields 
become increasingly cumbersome. For our purposes it is enough to present the expansion up to 
order O (e^), which is explicitly written in equations (A.l) of Appendix A.lp The companion 
expansion for the frequency is given by 



15 
28 



22456447 
35562240 



+ O (6«) 



(3.4) 



Having at our disposal the field expansion (A.l ) for {/, g, h, fi} we can now easily study their 
asymptotic fall-off and identify the relevant constants Cf,Ch,CQ introduced in (2.9). When 



plugged into (2.14) we find that the mass E and angular momentum J of the boson star are 



given, respectively, by 



E 



4 



5 2 77951 
6 ^ ^ 127008 



J 



TT 



1 2 83621 
12 ^ ^ 1270080 



+ O (6^) 



(3.5) 



Perturbative rotating boson stars are then a one-parameter family of solutions parametrized 
by the scalar amplitude. In particular, for small E, J, their mass is uniquely determined as a 
function of their angular momentum. (We will see in the next section that this does not remain 
true for larger E, J.) 



3.2 Hairy black holes 

Hairy black holes are a two-parameter family of solutions. In the perturbative regime, they can 
be parametrized by the asymptotic scalar amplitude e and, in addition, by the horizon radius 
r+. Consequently, to construct these solutions we need to do a double expansion in powers 
of e and also in powers of r+. Plugging this double expansion into the equations of motion 
we immediately find that, contrary to what happens in the boson star case, we cannot solve 
analytically the coupled system of perturbative equations everywhere. This occurs already 
at the linear level and is not a surprise since we already encounter a similar situation in the 
familiar studies of, e.g. the Klein-Gordon equation in the Kerr black hole. The traditional 
way to proceed is to use a matched asymptotic expansion procedure that applies quite well to 
the problem at hand, both at linear and non-linear level. In short, the main idea behind this 
approach is the following. We divide the exterior spacetime of our black hole into two regions; 
a near-region where r+ < r ^ £ and a far-region where r ^ r_|_. In each of these regions, some 
of the terms in the equations of motion make a sub-dominant contribution that is effectively 
discarded in a lower order expansion analysis. This yields a coupled system of perturbative 

^We should however note that when we compare our perturbative results of this section with the numerical 
results we will be using an expansion that encodes information up to order O (e^) . 
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equations that have an analytical solution in the spacetime region where the approximation 
is valid. The next important step is to restrict our attention to small black hole^ that have 
r+/£ ^ 1 (and e <^ 1). In this regime, the far and near regions have an overlaping zone, 
r+ ^ r <^ £, where the far and near region solutions are simultaneously valid. In this matching 
region, we can then match/relate the set of independent parameters that are generated in each 
of the two regions. Our matched asymptotic analysis of the rotating hairy black holes is similar 
in spirit to the analysis done in |T7] to construct static charged hairy black holes. 

We are ready to apply the matched asymptotic expansion to the construction of our hairy 
black holes. The first observation is that the expansion of the scalar frequency u must be 



independent of the region we look at since it defines the Killing vector (2.4) of the full solution. 
Its double expansion is given by 

n ^ 2i 

= X^^^2i(r+)e^^ uj2j{r+) = ^UJ2j,2i (^y) ' = UJ ] (3.6) 

jr=0 i=0 

i.e. at each order of the e-expansion we do a second expansion in r^/i. The last relation in 



(3.6) guarantees that the only Killing vector of the hairy black hole is also the null generator of 



its Killing horizon, and it is also required by regularity of the scalar field at the horizon; recall 



discussion associated with (2.12). 

The procedure now proceeds with the aforementioned split of the spacetime into the far and 
near regions. Describing the full details of this double expansion and matching procedure is in 
practice unfeasible due to the large expressions involved and intermediate steps required before 
matching/fixing the constants of the problem. In the next three subsections we will describe 
in detail the formalism and all the ingredients required to do the far-region, near-region and 
matching analysis. However we will not do the explicit analysis except in a small number of cases 
just to illustrate the computational task. The explicit final results for the expansions of all the 



fields up to an appropriate order is then presented in Appendix ( A.2 ). These expansions contain 



all the information we need to compute the energy, angular momentum and thermodynamic 



quantities in subsection 3.2.4 



Like the boson star, we will focus on only the ground state hairy black holes that have the 
lowest radial excitation of the scalar field. They form a two-parameter family that corresponds 
to the lowest possible energy to angular momentum ratio. The k^^-excited hairy black hole 
family of solutions could also be constructed following similar steps. 

3.2.1 Far-region expansion 

The far-region is defined by r 3> r+. In this region, a small hairy black hole is a small 
perturbation in r+/i and e around global AdS. 



To construct the hairy black holes, in addition to the expansion (3.6), we also do a similar 
double expansion in the gravitational fields that we denote collectively as F°^*, i.e. = 

and in the doublet scalar n°"*: 



n p 



j=0 i=0 



2i 



coordinate invariant definition of "small" black hole requires the black hole to have small energy and 
angular momentum. In the end of our analysis we will find that this definition is equivalent to have black holes 
with small r^/£ and e. 
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n 



out I 



r, e, r. 



2j+l 



Yfout 



(r, 



out / 
2j+l,2i\ 



j=0 



i ) 



2i 



(3.7) 



We use the superscript "out" to emphasize that these expansions are vahd only in the far-region. 



We want to find the gravitational and scalar fields that obey the equations of motion (2.6) and 



the asymptotic boundary conditions (2.9). The boundary conditions (2.11) at the horizon are 



not imposed since the far-region does not extend all the way down to this boundary of the full 
spacetime. 

Start by considering the leading order expansion in e, i.e. n = 0. Like in the boson star, 
this order describes the linear perturbation problem where we introduce a non-trivial massless 
doublet scalar in the spacetime background but this condensate does not back-react in the 
gravitational field. The difference with the boson star case is that now the background spacetime 
is not just global AdS but the {d = 5) MP-AdS black hole described in subsection 2.4 Since 
we know this solution exactly, the n = coefficients Fq 2* are already known up to any desired 
order 0{r^/i^). All we need to do is to insert the expansion (3.6) for uq into (2.21), that defines 
the mass-radius vm and rotation parameter a of the MP- AdS black hole in terms of r+ and 
flfj = uq, and replace the resultant expressions into the gravitational fields (2.19). A simple 



series expansion in r^/i gives the desired coefficients -Po'2i('^) terms of the constants u!o,2i 
(the latter are left undetermined at this stage). By construction, the associated F°"*(r, 0,r+) 



satisfies trivially the gravitational field equations (2.6a)-(2.6d) and the asymptotic boundary 
conditions (2.9) up to order 0{e,r^/i^). To find the coefficients 11°^*^ we now have to solve the 
scalar equation (2.6e) up to order C(e,r^J_/£^) subject to the asymptotic boundary condition in 



(2.9) for n. 



To illustrate the computational task let us discuss the details of the determination of the first 
two contributions n°'()*(r) and 111^2 {^)- The lowest order in the horizon radius expansion has no 
dependence on r+, i.e. it encodes no information about the existence of an horizon. Effectively 
we are determining the lowest normal mode in global AdS and the desired coefficients are 
precisely given by (3.3) after the relabelling {Fo,Ilo,uo} — )■ {Fq q*, cjq.o, n°"Q*}. The next-to- 



leading order is the first to encode a r+ dependence. We already described how we can get 
All we need to find is n°"2*- The requirement that the scalar equation is solved up to 



rpout 
^0,2 ■ 



C(e^,r^/£^) and the field decays as (2.9) demands that 



YfOut 
-■■^1,2 



24r3 (r2 + £2)^/2 
+6r^ (r^ + f) (2£wo,2 - 15) log 



- r(£wo,2 + 3) - 3r't{3iuo^2 + H) - 2r^t{Aiuo,2 + 45) 



r2 + ^2 



(3.8) 



Note that the asymptotic boundary condition was used to fix two integration constants that 
would give the more general but unphysical solution. At this stage the constant uo,2 is left 
undetermined. We will fix it in subsection 3.2.3 We leave a discussion of the linear analysis of 
the higher order terms in r+/£, as well as the discussion of the non-linear terms (!?(e"') (n > 2) 
to subsection 13.2.31 

At this point we still need to justify why the far-region analysis is valid only in the coordinate 
range r ^ r+. In subsection 3.2.3 we will take the small radius expansion of 11°"* and we will 
find that it is singular as r — )■ 0, diverging with a power of — . This shows that the far- region 
analysis certainly breaks down at r ~ r_|_. Therefore we must take its regime of validity to be 
r 2> r+. 
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3.2.2 Near-region expansion 



The near-region is defined by tlie radial range r+ < r ^ £. In tliis region, the rotational 
field Q ~ of the MP-AdS/hairy black holes is negligible compared with the mass scale 

£/r+ 1 set by the horizon radius. The best way to see this is to introduce rescaled time and 
radial coordinates, r = t£/r+ and z = ri/r+, such that the near-region is described by 2; ~ £. 
In these rescaled coordinates the d = 5 MP-AdS black hole (2.18)-( 2.19[ ) (to which we want to 
add a rotating condensate) is written as 



f 



h 




(3.9) 



{mj, - 1) - £2) ' rimi - z^ {ri {mi - 1) - p) ■ 

This coordinate transformation is illuminating because the rotation field of the conformal met- 
ric, Qrxi, (X ^Q, has an explicit factor of r+/£ <C 1. In the near-region the rotation field is 
therefore weak and it is to be seen as a small perturbation around the static configuration. It 
makes also explicit that an expansion in the small rotation field is is an expansion in a power 
series of r+/£. These statements should continue to hold when a small condensate is added to 
the system, i.e. when we do an additional power series expansion in e. 



These arguments justify that the expansion (3.6) is to be accompanied by a similar double 

{f 



expansion in the gravitational fields 



) y 



l^m^ and in the doublet scalar 11*' 



n N 

r-{z,e,r+) = ^F-(z,r+)e2^ i^^,"(^,r+) = Y.^ih^(^) ( 

7=0 i=0 



2i 



(3.10) 



N 



T-rm 



(z,r+) = ^n! 



out I 
2j+l,2j+l' 



j=Q 



The superscript "in" indicates that these expansions are valid only in the near-region. The 



near-region fields must obey the equations of motion (2.6) and the horizon boundary conditions 



(2.11). The asymptotic region z — )■ cxd is not part of the near-region so we do not impose the 



asymptotic boundary conditions (2.9). 



The strategy for the next steps of the near-region analysis is now similar to the far-region 
one. At leading order in the e-expansion, we again have the linear problem where we want to 
find the scalar field that propagates in the MP-AdS background without back-reacting with it. 
The gravitational coefficients Fq^^ are those of the MP-AdS background (3.9): we just replace 
the frequency expansion (3.6) into (3.9) and take a series expansion in r^/i up to the desired 
order 0{r^/£P). The non-trivial task is to determine the scalar coefficients ^T2i+i which satisfy 



(2.6e) up to order 0{e,r^/£P) and the last boundary condition in (2.11) 



As an example of the procedure, we find here the leading order contribution in the r^/i 
expansion of the linear scalar field. Equation (2.6e) is satisfied up to order 0{e,r'^/i'^) by a 
linear combination, {2z'^ / — l)-|-Q"Qi/2 {2z^ / P' — 1), of a Legendre polynomial of the 

first kind P and of a Legendre polynomial of the second kind Q. The latter has a logarithmic 
divergence at the horizon z = I. The boundary condition (2.11) therefore requires that we 
eliminate this contribution by killing its amplitude: = 0. The physical solution is then 

n- = (2^7^' - 1) , (3.11) 
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where the integration constant Kl^ will be fixed only by the matching analysis of subsection 



3.2.3 This is a good moment to note the curious fact that the far-region expansion of the scalar 
field is an expansion in even powers of r+/£, see (3.7), while in the near-region 11*" has instead 
an expansion in odd powers of r^/i, see (3.10). The reason why the 11*" expansion starts at 
r_|_/£ instead of at r'^/i'^ is that 11 must have a smooth limit when we send the horizon radius 



to zero, in order for perturbation theory to be consistent. One has n^"j^| 
n*" = e [ir{^-^r^/i + 0{r\/i'^)\ is regular when r+ — )■ as requirecfj 



— and thus 



L , , V ,j ^ , . Ultimately, this occurs 

because, as we will find in the end of our analysis, the boson start describes the r_|_ — )■ limit 
of the hairy black hole and in the former case one has n|r^o 



0; see (A.l) or (2.10). 



The discussion of the linear analysis of the higher order terms in r+/£, and the discussion 
of the non-linear terms (9(e") (n > 2) is left to subsection 3.2.3 A question we still have to 
address is why the near-region analysis is valid only for r <^ £. In subsection 3.2.3 we will take 
the large radius expansion of 11*" and we will conclude that it diverges asymptotically as j. 
Consequently, the near-region analysis breaks down at r ~ ^, and we must take its regime of 
validity to be r <^ £. 



3.2.3 Matching region 

Small black holes have r+ <^ £ and the far and near region solutions intersect in the region 
r+ <^ r <^ £ where they are simultaneously valid. Therefore, this overlapping region provides a 
set of matching conditions that fix the set of integration constants and frequencies 0J2j.2i- After 
this process the solution is fully determined everywhere. The matching procedure should be 
done iteratively order by order in the r^/i expansion for a given order of the e-expansion. That 
is, at a given order 0{e^r'^/£'^) we determine the associated unknown constants and only then 
we move to the next order 0{e''r'^^ /i'^'^'^) to do its matching analysis. The set of matching 
conditions are generated by matching the large radius behavior of the near-region solution with 
the small radius expansion of the far-region solution. To get the former, we take the large-z 
Taylor expansion of (3.10) and then we return to the original radial coordinate r = zr^/i. To 



get the small radius expansion of the far-region we do the small-r Taylor series of (3.7). The 



two golden rules for the matching procedure are the standard ones: at each order we should 
match only terms that will not receive a further contribution at next order of the expansion; 
and the divergent terms in the large radius near-region expansion (that go as powers of r/i) 
and the divergent terms in the small radius far-region expansion (that go as powers of r^/r) 
must be eliminated or at most exactly matched in the other region. 

Once more this matching procedure is better explained with an example. Consider the 
linear problem, i.e. the order 0{e). Here, we want to explicitly determine the scalar field up to 
order (9(er^/£^). Take first the order 0{er^/i). Insert 11™]^ as given by (3.11) into (3.10); do a 
Taylor expansion in z and replace z = r i/r^. This yields the asymptotic near-region expansion 



that we call 11]"^,^^. To get the small-r far-region expansion 11 

Tout 



small; do the Taylor expansion of 



n = en°Q. Matching these two expansions fixes the integration constant i^^" introduced in 



(3.11): 



IT*" 

large 



TJOUt 

■■--^small 



+ 



ir 



+ 0( 



large 



Y\out 

small 



-"■11 



71 

4 



(3.12) 



^To be more clear, 11* 



n 



1,1 '+/ 



o{rl/e'^)] with 



both the field equations and the horizon boundary condition. But it would not have a smooth limit when 
r+ — >■ and thus we have to set KIq = 0. The expansion then starts with the linear term in r+/i. 
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Take now the order 0(er^/^^). There is no coefficient ir{^2 near-region contribution to 

O(r^) comes exclusively from n*j^^^: we repeat the steps used to get (3.12 ) but this time we keep 



also the contribution in Hl^rge ^^^^ ^o be matched by the small-r far-region expansion. 

Matching the 



The latter is generated by the Taylor expansion of 11 



two expansions determines ix'0,2 introduced in (3.8): 



TTin 

large 



Y\out 
■'■■'■small 



4 £r 



+ 



1+ 



e 1- 



48 



24 



) +0 



, -TTin 

' large 



Y^out 
small 



1^0,2 



that we can now also insert in the expressions for the gravitational coefficients 
fixing them totally. We could now proceed and do the analysis at order (9(eri 



(3.13) 

0,: 

and 0{ert 



rpm J pout 

^0,2i '^i^*^ -^0,22 



The final results, after the matching, are presented in Appendix (A. 2). 
Having done the analysis at linear order 0{e 



contribution C(e'^ 



we must next consider the non-linear (n > 2) 
Indeed, only at these non-linear orders the scalar field back-reacts 



in the metric giving the necessary information about the gravitational field of the hairy black 



holes. At each order {n,p) we have to solve the coupled system of equations (2.6) at order 
O (e", r!J./£^). At any even order e", the equation for the scalar field (2.6e) is trivially satisfied 



and the gravitational field equations (2.6a)-(2.6d) are the ones whose solution encodes the non- 
trivial information on the back-reaction of the scalar field on the gravitational field. On the 
other hand, the equations at any odd order are used to find the corrections to the scalar field. 
The analysis again starts with the far and near region construction of the solutions that are 
then matched in the overlaping region. The process is done order by order and it is similar to 
the few cases we illustrated at the linear level. The analysis is however cumbersome for several 
reasons: the lower order contributions source the higher order fields; the expressions are quite 
often large; and there are many fields involved. For this reason we do not give the technical 
details here. The explicit final results for the expansions of all the fields up to combined 
order (n + p) 



4 is then presented in Appendix (A. 2). The near-region solutions with the 



matching analysis already done is presented in Equation (A. 3), while the far-region fields are 



written in Equation (A. 2), again after having fixed the integration constants of the problem by 
the matching procedure. The companion expansion for the frequency will be given below in 
equation (3.14). In the presentation of these results, we take e ^ 1 and r 
that 0{e] 



<^ i and we assume 

0{r^/i). The latter assumption implies that terms with the same {n+p) contribute 



equally to the perturbative expansion, i.e. O (e°,r^/£'^) ^ O {t 



~ C(e^rO/£0). We 



have done the consistent perturbative expansion analysis up to {n + p) = 4, which is sufficient 
for our purposes. Extending the perturbative analysis to higher orders is certainly possible, 
although increasingly cumbersome. 



3.2.4 Energy, angular momentum, and thermodynamic quantities 

The matched asymptotic expansion analysis of the previous subsections constructed pertur- 
batively the scalar and gravitational fields of the small rotating hairy black holes. With this 
information we can now compute their energy, angular momentum, and thermodynamic quan- 
tities and study some of their main properties. 

The angular velocity of the black hole, which is also the frequency of the doublet scalar 
field, was determined side by side together with the field expansions in the matched asymptotic 
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expansion construction; see eg (3.13). It is given by, 



5-3 
4 22456447 

— e 



959 3 , /r+ 
1^-2^°H4Z 
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35562240 
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(3.14) 



The perturbative near-region solution for {f,g,h} is written in Equation (A. 3) of Appendix 



(A. 2). It then follows from (2.15) that the temperature and entropy of the rotating hairy black 



holes are, respectively, given by 



H 
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(3.15) 



Given the perturbative far-region solution (A. 2) for {f,g,h,il} we can read their asymptotic 



fall-off and identify the relevant constants Cf,Ch,Cn introduced in (2.9). Plugging these into 



(|2.14|) we find that the energy E and angular momentum J of the rotating hairy black hole are 

4 / 77951 



given, respectively, by 



E 



J 



TC 
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(3.16) 



It can be explicitly checked that relations (3.14)-(3.16) obey the first law of thermodynamics 



(2.17) up to fourth order. 



To discuss these results, start by noticing that the perturbative analysis constructs small 
black holes. Indeed r+/^ <^ 1 and e ^ 1 implies that E/i"^ <^ 1 and J/i^ ^ 1. One of the main 
properties of the hairy black holes is that they reduce to the boson stars constructed in section 
3.1 when the horizon radius is sent to zero. A priori there is no reason why this had to be 



the case. The boson stars may have been analogous to marginally bound states which would 
disappear when a small black hole is added. However, we indeed find that solutions exist and 



the scalar field frequency (3.14) and the conserved quantities (3.16) reduce to the frequency 
(3.4) and energy and angular momentum (3.5) of the rotating boson stars in the limit — )■ 0. 



Naturally, the black hole entropy also vanishes in this limit, while its temperature diverges (a 
similar situation — see (2.22) — occurs when we send r+ — > in the MP-AdS black hole). 
Another special curve is the one-parameter sub-family of hairy black holes with e = 0. In this 
limit the scalar condensate vanishes: this curve describes the merger of the hairy black holes 
with the MP-AdS black holes. This merger line describes non-extreme MP-AdS black holes 



that are at the onset of the superradiant instability, as will be discussed later in section 4.1.1 



In sum, in a phase diagram oi E/C."^ vs J/£^ the two parameter family of hairy black holes is 
bounded below by the curve describing the one-parameter family of boson stars, and bounded 
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above by the merger curve with the MP-AdS black holes. We will postpone further discussion of 
the properties of the system to section |4] where we will construct numerically the most general 
family of rotating boson stars and hairy black holes whose energy and angular momentum 
are not restricted to be (very) small. We will find that the phase diagram develops a very 
interesting structure as the energy and angular momentum of the boson stars and black holes 
grow to larger values. 



4 General results 



The perturbative construction of the small rotating boson stars and hairy black holes of the 
previous section has the advantage of being an analytic analysis. However, as we will find 
latter in this section, there are interesting properties of the system that can only be revealed by 



studying solutions of (2.6) in regimes where the energy and angular momentum are not small, 
i. e. where the analytical results are not valid. To probe this region of the parameter space we 



thus need to perform a frontal numerical approach in solving the equations of motion (2.6). 



This will be the aim of this section. As a check of our numerical work, for small energy and 
angular momentum we will find agreement with the analytical results found in section [3} for 
both the boson stars and black holes. 



(2.8). 



Before proceeding, we recall that the metric ansatz (2.2) has the residual gauge freedom 
We will use this freedom to optimize the numerical construction of our solutions. In 
this computational task, we either assume A to be such that f2(oo) = (the physical gauge, as 



explained in subsection 2.2) or u = 0. Note that the latter necessarily implies that the bound- 
ary metric is rotating, which corresponds to a rather unphysical situation. Any calculations 
developed in this gauge, should (and will) therefore be recast in a gauge where f2(cxD) = 0, 
which in turn implies u ^ 0. It turns out that these two gauges are not equivalent numerically, 
in the sense that convergence is achieved at different rates in different regions of the parameter 
space whether we use one or the other. 

In constructing both the rotating boson stars and the rotating hairy black holes we use 
a numerical method which relies on a standard relaxation procedure, combined with spectral 
methods defined on a Chebyshev grid. In order to solve the coupled system of differential 



equations (2.6), we first discretize Eqs. (2.6a 2.6e) and use the constraint equations, Eqs. (2.6f 



2.6g), to gather information about the appropriate boundary conditions at asymptotic infinity 
and either horizon (black holes) or origin (boson stars). Once the information about the 



constraints is imposed on the boundary conditions (2.9), (2.10) and (2.11), the flow of the 
equations of motion into the bulk ensures that the constraints are satisfied everywhere, see 



Eqs. (2.7) 



A brief but incomplete outline of this section follows. We start by presenting the growth 
rates of the superradiant instability for several values of the azimuthal quantum number m. 
Not surprisingly, we will find that the onset of this instability signals, in a phase diagram, a 
merger or bifurcation line that connects the MP-AdS black holes with the rotating hairy black 
holes (one merger line for each m). This conclusion is firmly established once we compare the 
m = 1 threshold instability curve with the hairy black hole curve we obtain independently 
when we send the condensate to zero in our hairy black hole solutions. We find that the most 
unstable mode is the one with m = 1. As discussed previously, the cohomogeneity-1 ansatz 
(2.2), (2.3) is appropriate to study the non-linear solutions that bifurcate from the m = 1 
merger line. We thus use this ansatz to construct the rotating boson stars and rotating hairy 
black holes, this time numerically and for any values of the energy and angular momentum. 
This section will culminate with the presentation of a full phase diagram of solutions where we 
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gather the information on the m = 1 hairy solutions (both numerical and perturbative results) 
and on the MP-AdS black hole. The hairy solutions we present here are associated with the 
scalar superradiant instability. However, MP-AdS black holes are also superradiant unstable 
against gravitational perturbations. We will also study the growth rate of this gravitational 
instability (the onset of the instability was first studied in [1]), and find that the properties of 
this instability are similar to those of the scalar instability. 



4.1 Linear expectations 



Here, we intend to provide three distinct properties about rotating hairy black holes solutions 
that branch off the MP-AdS black holes. First, we determine the locus in the moduli space of 
solutions where the hairy black holes merge with the MP-AdS family. Next, we investigate the 
growth rate of the superradiant instability associated with the scalar field H, in the background 
of a MP-AdS black hole. Finally, we will find that the properties of this scalar instability are 
similar to those of the purely gravitational superradiant instability. 

The equation we want to solve is that of a massless scalar field doublet, H, in the background 



of the MP-AdS solution, described in section 2.4 

V^II = 0. (4.1) 
Because we would like to compare the growth rates associated with different values of the 



azimuthal quantum number m, we need to generalize (2.3). Even though we did not manage 



to extend the inclusion of m > 2 modes in the construction of the non-linear solutions, at the 
linear level this generalization is easily achieved: 



II = H(r) e 



-i LU t+i m tp 
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(4.2) 



Inserting this expression in (4.1), we obtain the following ordinary linear differential equationj^ 



+ 



h{ijj — mVlY 2\m\ 



f 



m 
r'^h 



H = 0. 



(4.3) 



In order to solve this equation we need to discuss first the boundary conditions to be 
imposed. Naturally, we are interested in ingoing boundary conditions on the future event 
horizon, which means that close to the horizon r = r+, H(r) has the behavior: 



U(r) 



Exp 



2 + 2rl (1 



r\ + 



^log 1 



ao + C 1 



_ !+ 

(4.4) 



where is an arbitrary constant, and we defined 

= a; — m£Lu . (4.5) 
On the other hand, at infinity there is one normalizable mode only, for which H decays as 

W)\.,^^^0{t-'), (4.6) 

like the non-linear solution. 

^This expression matches expression (33) in for /i = fc = 0, = (t = 1 and — r'l'^\fhX\. 
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4.1.1 Merger 

As we mentioned above, the rotating hairy black holes will branch-off the MP-AdS solution 
for specific values of the MP-AdS angular velocity Qh and horizon size r_|_. The strategy to 
determine this zero-mode is simple. We know that u = mVtii (w^ = 0) for all hairy black holes: 



see discussion associated with (2.12) for the particular case m = 1; the result for m 7^ 1 follows 
from a similar reasoning. In particular, this statment is valid for the merger line (where the 
condensate vanishes) , so we can choose values for both m and r_|_ /£ and ask what is the value 



of VLh^ for which the solution exists. Unfortunately, we did not manage to write (4.3) as an 
eigenvalue problem in VLhI, so we resort to a shooting procedure using a standard fourth order 
Runge-Kutta method. This can be easily implemented in the following way: 1) we choose 
a pair (r+ / £, m) and determine the solution near the horizon and in the asymptotic infinity 
region using the Frobenius method; 2) then we integrate from both boundaries to a point in 
the middle of the integration domain where we demand the two solutions and their derivative 
to agree; 3) this, in turn, will only occur for a specific value of ^2^^, which we determine using 
Newton's method. 

Each value of m labels a different merger line. We are primarily interested in the m = 1 
mode because we will find that these are the most unstable, and the associated black holes are 
those we construct at the full non-linear level. This is the first step to determine the full phase- 
space of solutions, since it provides a rough estimate of which energy and angular momentum 
scales we need to sweep. The results are illustrated in Fig. [1} We see that Qh^ starts being five, 
for very small black holes, in agreement with the analytical results for the radial mode k = 
of the previous section, and slowly decreases to Qh^ ~ 4.8652, for extremal MP-AdS. These 
results are in very good agreement with the analysis performed in [1]. The m = 1 merger line, 
here plotted in the Qh^) plane, can equivalently be expressed in terms of the energy and 



angular momentum of the MP-AdS solution using Eqs. (2.23). This is the representation we 
will use in sections 14.31 and 14.41 

The other two merger lines in Fig. [T] correspond to the merger lines of the m = 2 and m = 3 
modes. They share most of the characteristics of the m = 1 mode, except that the critical 
angular velocity above which small black holes are unstable changes. For m = 2 this happens 
around = 3 while for m = 3, ^Ih^ = 7/3, in agreement with Also, the higher the value 
of m, the larger the unstable MP-AdS black holes can be, because the unstable region extends 
all the way up to extremality. The limiting value of fin as r+ — )■ is directly related to the 
modes of a massless scalar field in AdS^. The lowest radial mode with angular dependence m 
{i.e. having eigenvalue m(m + 2) for the Laplacian on S^) has frequency ui = 4 + m [18]. It is 
thus invariant under the Killing vector dt + {u/m)d^. If we want to add a small black hole, it 
must have angular velocity Qh^ = w£/m = 1 -|- 4/m. This reproduces the values 5, 3, 7/3 that 
we find for m = 1,2, 3. The point is simply that the mode of the scalar field in AdS^ is the 
marginal mode we are looking for when the black hole is very small. 

When Qh^ < 1 the MP AdS black holes do not have a superradiant instability for any m. 
This is because the Killing field dt + ^nd^ remains timelike everywhere outside the horizon. 
Using this symmetry one can show that there is no superradiance for any field satisfying the 
dominant energy condition [TT] . 

For values of the angular velocity above the merger lines, we expect the MP-AdS solution 
to be unstable. This is precisely what we are going to confirm in the next subsection. 
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Figure 1: Plot of angular velocity versus horizon size. MP-AdS black holes exist in the shaded 
regions and its boundary corresponds to extremal black holes. Curves of f2//(r+, m) correspond- 
ing to the merger lines are displayed for different values of m. Above each line, black holes are 
unstable to perturbations with that m. Below VLh^ = 1, black holes are stable. The inset plot 
zooms in the region where the m = 1 merger exists. 



4.1.2 Growth rates 



At this point, we know exactly where the rotating hairy black holes branch-off the MP-AdS 
solution. Now, we want to quantify how unstable the MP-AdS black holes are. For this we 
need to provide the growth rates of the superradiant instability which afflicts the MP-AdS 



black holes. Namely, we need to go further in solving (4.3) and consider solutions for which 



u 7^ mVlH- In order to work with analytic functions only, we define the auxilary function q as 
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(4.7) 

where we have implicitly performed the change of variables y = 1 — r^/r^. Note that y is a 
compact variable taking values in the unit interval. By construction, q{y) is an analytic function 
both at the future event horizon and at asymptotic infinity. Substituting this expression into 



(4.3), gives the following quadratic Sturm-Liouville problem 

Log - (^*jC-iq - ujlC2q = 0, 



(4i 



where each Ci is a second order differential operator in y, independent of cj*. The procedure 
is now clear: we first transform q{y) into a vector, say q, where each of its entries is given by 
q{y) evaluated at the Chebyshev points yj. In this discretization scheme, the operators £j are 
matrices, say Lj, that act on q. Consequently, equation (4.8) reduces to a generahzed eigenvalue 
problem in w*, 

0, (4.9) 
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where I is the identity matrix. This linear system can be easily solved by the in-built Mathemat- 
ica routine Eigensystem. In order to generate the growth rates, we can choose values for r+/£. 
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Figure 2: The real (left panel) and imaginary (right panel) parts of the perturbation frequency 
as a function of angular velocity, for a constant value of the MP-AdS horizon size, r_|_/£ = 0.1. 
Different curves correspond to different values of m. Note that the plot on the right uses a 
logarithmic scale. 



VLh^ and m and then compute using (4.9). Note also that there are no reality conditions 
on o;^. Therefore, shooting methods would certainly face serious challenges in solving (4.3) for 
generic values of r+/£, VLh^^ and m, because cj* is generically complex. 

The results we obtained were generated along a line of constant MP-AdS horizon size, 
r+/£ = 0.1, and are illustrated in Fig. |2} As a check, we note that the value of VLh(- for 



which Ivniujt) = is precisely the zero-mode described in the Merger subsection 4.1.1, The 



concordance between the onset critical angular velocity obtained using the shooting method of 
the previous subsection and the spectral methods of this subsection is excellent; in particular, 
they agree with each other within a 0.001% error. 

The first striking property about the growth rates in Fig. [2] is that they are all small, 
leading to a secular type instability. This is reminiscent of what happens in the case of a 
massive scalar field in the background of the asymptotically flat Kerr black hole [in]-[23]- Also, 
as m increases, the growth rates become increasingly smaller. In fact, each unit of m seems 
to change the growth rate by two orders of magnitude. These results indicate that the m = 1 
mode is the most unstable mode, when coexisting with the remaining ones, and that modes 
with higher m will become increasingly less unstable, in the sense that their maximum growth 
rates become increasingly smaller. The fact that the real part of cj^ is negative just means that 
< Re{u) < mVlHi which is a well know property of the superradiant instability |24j-|27]. 



4.1.3 Linear gravitational superradiance 

Before proceeding to the numerical construction of the rotating boson stars, we would like 
to compare the scalar instability described in the previous subsections, with the superradiant 
instability in a purely gravitational setup. It is well know that Kerr-AdS black holes, and 
their natural extension to higher dimensions p2j, exhibit the superradiance instability [TT]. 

In this subsection, we will argue that much can be understood about this instability 
by inspecting our (more tractable) cohomegeneity-1 system. In particular, we will show that 
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Figure 3: The real (left panel) and imaginary (right panel) parts of the perturbation frequency 
as a function of angular velocity, for the purely gravitational superradiance instability, at fixed 
MP-AdS horizon size r+/£ = 0.1. The mode shown has ii = m = 3. 



the purely gravitational instability has a spectrum which is very similar to the one in Fig. [2j 

To study this problem, we would need to perturb the five-dimensional MP-AdS line element 
described in section |2.4[ This, however, turns out to lead to a complicated set of coupled 

which will be studied elsewhere. Here we take a simpler 
The reason why it is so challenging to perturb the five- 



ordinary differential equations [33] 
route which suffices for our purposes, 
dimensional MP-AdS metric is related with the fact that no tensor harmonics exist on the 
CPi = S"^ base space, which in turn implies that we have to deal with either vector-type or 
scalar-type perturbations. In order to avoid this conundrum, we change gears and consider the 
seven-dimensional counterpart of this geometry, which was studied in The d = 7 equal 
angular momenta MP-AdS is still cohomogeneity-1 (see subsection 2.4[ ) and its base space is 
CP2, where tensor harmonics do exist. This is still a reasonable comparison because we do not 
expect the physics involved here to depend on the number of dimensions. The equation for the 
tensor-mode was derived in |^, and can be written as: 



J_ 



f 



where 

V{r) 



f 



.5/2^3/4 



Vh 



5/2V 



h 



(4.10a) 



r'^h 



£i(£i +4) -m^ 



(4.10b) 



and / and g are given by (2.19) with N = 2. Here, \l/ is the master variable and ii is a positive 
integer that can take certain values that depend on m. We will choose the values ii = m = 3, 
which were studied in [3]. Also, for the sake of concreteness, we will only consider a line of 
constant MP-AdS horizon size, namely r^/i = 0.1. The authors of |1] studied the merger line 
associated with the gravitational superradiance. Here, we will go further and compute also the 
instability growth rate. 

We present the results in Fig. [3} where we can see that the purely gravitational setup shares 
most of the features of the scalar instability described in the previous subsections. There is, 
however, a quantitative difference, namely, the instability growth rate seems to be even smaller 
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than the one associated with the doublet complex scalar field, at least in the tensor-type sector 
of the perturbations. It is not clear whether this difference arises because we are comparing a 
five-dimensional result to a seven-dimensional result or because tensor-type perturbations are 
generically smaller. Another important difference, which we will address next, is that for the 
complex scalar field doublet we can find the endpoint of the instability in the m = 1 sector. 
The gravitational case seems to require solving PDE's which are not of the standard elliptic 
type (since the Killing field becomes spacelike at large radius). A full time-dependent evolution 
may be needed for any type of gravitational superradiant perturbation. 



4.2 Rotating boson stars 

We now turn to the numerical construction of boson star solutions to (2.6). Recall that rotating 
boson stars form a one-parameter family of solutions. In the perturbative construction of section 



3.1 the scalar amplitude e was an appropriate quantity to parametrize the solution. This is no 
longer the case when we extend to larger values of the energy and angular momentum, since e 
no longer defines uniquely the solution. A better choice is the energy density at the center of 



the boson star. For the ansatz (2.3), the energy momentum density at the centre of the boson 
star turns out to be proportional to n'(0)^, and as such we will parametrize our solutions by 

n'(o). 

Before proceeding, a few words about the numerical construction of these solutions are in 
order. Since we are using relaxation methods on a Chebyshev grid we need to guarantee that 
the functions we solve for are analytic in the domain of integration. In order to achieve this, 
we define the following auxiliary functions {qf^qhAnAg^lu}'- 



„2 



air) = 1 + n(r) = jt^I^ U 11) 

y^^^ ^+(^2 + £2)2' (^2+£2)5/2- ^^"^^^ 



For numerical purposes it is also easier to work in a gauge where the line element (2.2) is 
rotating at infinity, and for which u = 0. In the end we will change gauge, and set fi(oo) = 0, 
thus generating a non- vanishing value for u. Finally, in using spectral methods it is convenient 
to compactify the integration domain, which can be readily achieved by introducing the new 
radial variable, 



whith ?/ G (0, 1). At y = 0, regularity fixes the boundary conditions (2.10), and at y = 1 
we demand the spacetime to be asymptotically global AdS, i.e. that it obeys (2.9). These 
boundary conditions imply the non-trivial relations: 

qfiO) = 0, g,(0) = , ^^(0) - 2g,(0) - q',{0)[l + g,(0)] - ^[1 + q^H = , 

fo(0) - |gn(0) = 0, gn(0) - go = , (4.13a) 



and 



g}(l) + qf{l) + g,(l) = 0, g,(l) + g,(l) = 0, ^^(1) + 2g,(l) = , 

q'^il) = 0, q'M + ]^[25 - gn(l)1gn(l) = , (4.13b) 
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where go = ^n'(O) should be regarded as a number that we vary as we move along the boson 
star family of solutions. 



In Figs. 4(a) 4(d) we show how the dimensionless energy, angular momentum, frequency 
and asymptotic scalar field amplitude change as we increase the energy density. The solid lines 
represent the exact numerical solution, whereas the dashed lines show the analytical estimates 



of section |37L} The physics involved here is in many respects similar to the fiat space counterpart 
[6]. We see that there is a sign of a damped oscillatory behavior, similar to the charged solitonic 
solutions in AdS^ x 5*5 found in |3l]. This behavior occurs around a non-zero critical value of 
the energy, angular momentum, scalar frequency and scalar amplitude, respectively. We will 
label these values by {Ec, Jc, uJc, ^c}, respectively. The analytic calculations developed in section 
3.1 are in very good agreement with our numerical data in their regime of validity, i.e. for small 



enough n'(0). Note however that the perturbative analysis does not capture the interesting 
phenomenon of the damped oscillations. 

There is a simple consequence of the damped oscillatory behavior. Consider for example 
a plot of the condensate strength e as a function of the scalar field frequency u (we could 
alternatively choose other quantities, but the general picture would be unchanged). Due to the 
damped oscillations of Figs. 4(c) |4(d) , we can view e and to as coordinates of a two dimensional 
damped harmonic oscillator. Typical motion in the e — u plane is a spiral and that is indeed 



what we find in Fig. 4(e 



It is now clear that the line of boson stars ends at a specific energy and angular momentum, 
namely at Ec and Jc, respectively. As we observe in the top two panels of Fig. [4], we would need 
an infinite energy density or, equivalently, an infinite n'(0) to reach the critical values Ec and 
Jc- A possible detection of this singular behavior can be readily observed by computing the 
Kretschmann scalar Kg at the center of our geometry. This quantity, in terms of the spectral 



functions defined in (4.11), is simply given by 
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^X, ^ i'R''""'Rabcdl^, = 40 - 32 g2 + y go' + 96 q'M- 



(4.14) 



Fig. |4(f)| shows how the Kretschmann scalar evolves with tui. In particular, it seems to diverge 
at Uci, in accordance with our expectations. 



We expect an infinite number of damped oscillations in the plots of Figs. 4(a) 4(d), as the 
system evolves along the curves towards the critical values {Ec, Jc,uJc,^c}- As a consequence, 
the associated spirals are predicted to continue ad aeternum in Fig. 4(e) In particular, this 
means that for a line of constant u = Uc, we expect to find an infinite number of solutions with 
different values of J, say. The reason why we do not see the spirals continuing further inwards 
in our plots {i.e. why we do not observe the other arms of the spiral) can be well understood 
by inspecting what is happening to the curvature invariant Kg. Indeed, Fig. 4(f) clearly shows 



that the curvature is quickly diverging as we move inward along the spiral. Unfortunately, this 
makes it considerably harder to find numerically the solutions. 



As a check on our numerics, we have verified that our solutions satisfy the first law (2.16) 



with only 0.005% error. This provides a measure of the error in determining these solutions. 



4.3 Rotating hairy black holes 



The perturbative construction of section |3.2| necessarily explores only the small mass and an- 
gular momentum region of the parameter space of hairy black holes. In this section, we want 
to go beyond this perturbative construction, and construct hairy rotating black hole solutions 
with a single Killing vector field whose energy and angular momentum are not restricted to 
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Figure 4: Some properties of rotating boson stars. The solid curves represent the exact solution, 



whereas the dashed lines show the analytical estimates of section 3.1 
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be small. The construction procedure will be very similar to the one followed in the previous 
subsection for the boson stars, except that the black hole solutions will now depend on an extra 
parameter, namely the black hole horizon size r+. It is tempting to parameterize our solutions 
by r+ and Qh- Indeed these two quantities uniquely label the MP AdS black holes. However, 
we have already seen that for the boson stars u is not single valued, and we know that the 
hairy black holes must have Qh = ^- This strongly suggests that the hairy black holes will not 
be uniquely described by r+ and Qh, and this is indeed the case. 

It turns out that there is a simple geometric way to uniquely label the hairy black holes. 
These five-dimensional black holes have squashed horizons, which can be viewed as 
bundles over S"^. A natural parametrization of our solutions is just given in terms of the size 
of each of these spheres, which we will label by ri and r2, respectively. In terms of the metric 



ansatz (2.2), these are given by 

ri = r+ A//i(r+) , r2 = r+ . (4-15) 

We will explore this two-dimensional parameter space by fixing various values of r2 and slowly 
increasing ri. 

Numerically, this is difficult to implement directly. To input ri, we need to work in a gauge 
where a; = 0. However, we then find that "half way" through our computations, the code 
develops a numerical instability. This can be cured by changing to a second (different) code. 
In this new code, we work in a gauge where f2(oo) = and input u directly. This makes 
the construction of these black holes challenging, because each line of constant r2/i needs at 
least two independent codes to be generated. In addition to this, all the black holes we found 
are rather small, namely r+/£ < 0.14, which means that the curvature scales involved in the 
construction of these solutions are large. For these reasons, we find appropriate to use a rather 
dense Chebyshev grid with no less than two hundred points. 

For numerical reasons, it is convenient to introduce the new functions {qf,qh,qg,qn,(ln} 



which are adapted to our boundary conditions (2.9) and (2.11) 



m = ( S + iVl - 41 Qf{r) , h{r) = q,{r) , g{r) = qg{r) , H(r) 



^2/'^/VW5 '"VW J , tlV / igv J , (-^2 _^ £2^5/2 ' 

(4.16) 



1 / \ 

^^^\r) = - 1^1 - -± J gf^(r) , if gauge u = 0; Q'^'^\r) = ^ ' S^^S^ ^(°°) = ^■ 

We emphasize that f2(r) will either be given by fi^^^(r) or ^l^'^\r), whether we use a gauge 
where a; = or f2(cxD) = 0, respectively. Moreover, we use the radial variable ?/ = 1 — r^/r^, 
which takes values in the unit interval, the horizon being at |/ = and asymptotic infinity 
at y = 1. For y = 1, we demand the spacetime to be asymptotically AdS, and at ?/ = we 



require regularity. The system then obeys the boundary conditions (2.9) and (2.11). These, 
in turn, demand certain non-trivial relations among the functions {<?/, Q'g, Q'n, Q'q} and their 
derivatives evaluated at y = 1 or y = 0, four of which are similar in many respects to the boson 



star boundary condition relations (4.13). However, for the black holes, these four boundary 



conditions cannot be expressed in such a compact form as (4.13), and as such shall not be 
presented here. The remaining boundary condition imposes ^{r^) = u and ensures regularity 
at the black hole horizon. 

Our findings are compiled in Figs. |5| Each panel contains a set of lines that label different 
values of r2/i. From left to right in each panel, one has r2/i = 0.05, 0.07, 0.09, 0.11, 0.13. Let 
us start with the two top panels, where we observe the variation of the energy E and angular 
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momentum J with ri / ^. For small enough r2 / we observe an oscillatory behaviour that seems 
to stop abruptly when we reach extremality; see Fig. 5(c)[ Furthermore, as r2/£ decreases 
to zero, the number of oscillations before reaching extremality seems to increase. This seems 
to make sense, since we expect an infinite number of oscillations for the boson stars, i.e. as 
r2 = r+ — 0. In Fig. 5(d) we show how the black hole angular frequency Vt^ = u varies with 
ri/i and, as anticipated, the frequency starts precisely at the m = 1 zero-mode determined in 
the Merger subsection 4.1.1[ see the m = 1 curve in Fig. [1} The relat ive error in determining 
the MP-AdS onset angular velocity using the linear results obtained in subsection |4.1.1[ and 



the full non-linear results of this section (when the hairy black holes are close to the MP-AdS 
solution) is smaller than 10~^%. So, the onset of the m = 1 superradiant instability signals, in 
a phase diagram, a merger or bifurcation line that connects the MP-AdS black holes with the 
m = 1 rotating hairy black holes. Also quite interestingly, the black hole angular velocity is 



bounded below by the smallest frequency of the boson star: compare Fig. 4(c) with Fig. 5(d) 
In particular, we find no hairy black hole with Qh"^ < 1- 

We can also study the nature of the extremal hairy black hole solutions, in particular whether 
they are singular. In order to test this, we first studied the Kretschmann scalar evaluated at 
the horizon as a function of ri/i and r2/i. However, we soon realised that the Kretschmann 
scalar remains finite as we approach extremality. This means that if the solution is singular, 
it cannot be detected by this curvature invariant. Fortunately, in [35j, the authors proposed a 
new mechanism to test whether a black hole is singular. In particular, they observed that some 
large static black holes have all curvature invariants finite at the horizon, yet for such black 
holes, any object which falls in experiences enormous tidal forces outside the horizon. In order 
to check whether the hairy black holes are afflicted by the same problem, we need to calculate 
the geodesic motion of a timelike particle around our solutions. We label its geodesic tangent 
vector by X"'. This turns out to be a rather simple exercise, because the form of the metric we 
are currently studying has enough symmetry to reduce this problem to a set of quadratures. In 
particular, we consider a radially infalling particle in the equatorial plane and initially at rest 
at asymptotic infinity, for which 




(4.17) 



In Fig. 5(e) , we plot 



=r+ as a function of ri/i, and we see that it diverges precisely 
as we approach extremality, signalling a singular black hole solution, i.e. rotating extremal 
hairy black holes are singular. 

Finally, we point out that the hairy black holes solutions obey the first law of thermody- 



namics (2.17) up to 0.005% error. Like in the boson stars, we view this as a very good check 



on our numerics, and as providing a measure of the error in determining these solutions. 



4.4 Complete phase diagram 

In this section, we will sew all of the information gathered in the past two subsections, and 
construct the phase diagram for these rotating hairy solutions, i.e. the phase diagram for MP- 
AdS black holes, m = 1 rotating hairy black holes and rotating boson stars. As we shall see, 
the phase diagram is rather intricate, but the overall picture is quite simple. 

We will present the results in the microcanonical ensemble. In Figs. [6} the dashed blue line 
represents the MP-AdS extremal line, above which non-extremal regular MP-AdS exist. In 
both plots, the solid purple line represents boson stars. 
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Figure 5: Some properties of rotating hairy black holes. Each line corresponds to a different 
value of r2/£ = r+/^. From left to right in each panel, one has r2/£ = 0.05, 0.07, 0.09, 0.11, 0.13. 
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In Fig. 6(a) we show how the boson star energy varies with its angular momentum. One 
can understand the shape of this hne as follows. As we have seen in section 4^, the energy E 
and angular momentum J of the boson stars oscillate as a function of the energy density at the 
centre of the star, n'(0). However, unlike Fig. 4(e)[ the first law ensures that E and J have 
extrema at the same values of n'(0), which we label by n'(0)j. This means that instead of a 
spiral, the motion in the E, J plane has cusps. Let us define Ei{J) to be the energy of the star 
as a function of J for the branch where n'(0) e (n'(O)j, n'(0)i+i), where we identify n'(0)o ^ 0. 
We can see up to two cusps (and thus three branches of boson stars), one of which is plotted as 
an inset plot. Furthermore, we find that some boson stars coexist with MP-AdS black holes. 
This happens when the solid purple line is above the dashed blue line. 

The hairy black holes lie in a thin band above the boson star curve. To make this clearer, 
it is convenient to plot the energy relative to the lower line of the boson stars, i.e. AE = 



E{J) — Eq{J). This is shown in Fig. |6(b) The dashed blue line now returns to zero at the 



point where the boson star crosses the extremal MP AdS line in Fig. 6(a) The solid purple 



line on the far right is the second branch of boson stars. The green "vertical" curves represent 
hairy black holes with constant angular frequency ^h^-, whereas the dotted red curves are lines 
of constant r2 = r+. Recall that ri is increasing along the lines of constant r2- Thus, whenever 
these dotted curves cross we have hairy black holes with different horizon geometry but the 
same energy and angular momentum. In other words, we have non-uniqueness. Since there is 
also a MP-AdS solution, we have up to three coexisting black hole solutions. 

Close to the merger, the entropy of the MP-AdS solutions is smaller than that of the hairy 
black holes, being equal exactly at the merger line. This means that this transition is a second 
order phase transition. However, for sufficiently large angular momentum, the MP-AdS black 
holes also coexist with the hairy solutions, their entropy being now larger than that of the hairy 
black holes. Moreover, the transition is now first order, because these solutions never merge for 
this range of angular momentum. In sum, in a 3-dimensional plot of {S/£^, AE/i"^, J/i^} the 
hairy black hole family is a 2-dimensional surface bounded by the merger line and the boson 
star curve, if the angular momentum is smaller than the critical J where the red and blue 
curves meet in Fig. |6| This surface continues for larger values of J where it is now bounded by 
the boson star line and the extremal hairy black hole curve. This 2-dimensional surface never 
intersects with itself and has a sequence of (regular) "cusp lines" . 

We have provided evidence in section 4^ that the zero temperature limit of the hairy black 
holes is singular due to the appearence of large tidal forces. In order to estimate the shape 
of this singular curve in the phase diagram, we take the points of smallest temperature for 
each line of constant r2/i and Qh^- This line is represented as a dotted-dashed green curve 
in Fig. 6(b) The entropy along this curve monotonically decreases as we move away from the 
merger line. We conjecture that this line should join to the endpoint of the boson star, which 
has zero entropy, perhaps spiraling an infinite number of times before reaching it. 

To end this section, we will compare the numerical results with the analytical estimates 
of section 3A and 3^ In Fig. [7] we plot the analytical estimates and our exact numerical 



data. The agreement is excellent for small E and J. We might have expected our perturbative 
analysis to be valid for a larger range of energy and angular momentum, because the black 
holes are always small. However, we point out that the curvatures involved here are rather 
large {i.e. the coefficients of the expansion get quickly large as the order increases), and thus 
these expansions should have a small radius of convergence. This is the reason why we observe 
large discrepancies in the analytical estimate of the merger line for J/i^ > 0.001. 
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0.00 0.02 0.04 0.06 0.08 



(a) Energy of the boson star (solid line) as a function of its 
angular momentum. The inset plot shows {E2 — Ei)/P as a 
function of J/l^. {Ei,E2 are defined in the main text). 




0.00 0.02 0.04 0.06 0.08 

(b) Complete phase diagram for rotating hairy solutions. Ai? 
is the energy relative to the lower branch of boson stars. 

Figure 6: Phase diagrams of rotating hairy solutions. The dashed blue curve indicates the 
location of the extremal MP-AdS black holes, and regular MP-AdS black holes only exist above 
it. The solid purple curve represents the boson stars, the dotted red curves represent lines of 
constant r2/£ = Tj^/C. and the green "vertical" solid curves are lines of constant VLh^- Finally, 
the dotted-dashed green curve represents the best approximation to the singular extremal hairy 
black holes. 
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J/f" 

Figure 7: Phase diagram for small values of the energy and angular momentum. From top to 
bottom, analytical estimate for the merger line, exact merger line, extremal line of the MP-AdS 
black holes, analytical estimate for the bosons stars and finally the exact data for the bosons 
stars. 



5 Discussion 



By exploiting a clever ansatz for two (complex) scalar fields [6], we have constructed hairy 
black holes with only a single Killing vector. The ansatz corresponds to an m = 1 mode of 
the scalar fields. When the horizon area shrinks to zero, the solutions reduce to rotating boson 
stars. These solutions were constructed using two complementary approaches: a perturbative 
analytic approach for very small black holes, and numerical methods for larger ones. The phase 
diagram of these solutions was given in Fig. 6 and it was found that they exhibit nonuniqueness. 

We now make a few comments about the stability of our solutions. First consider the boson 
stars. We know that AdS is linearly stable, and because boson stars with small E and J are 
just perturbations of AdS, we expect them to also be hnearly stable. In fact, we expect the 
normal modes of these boson stars to be simple deformations of the AdS normal modes (3.2). 
For larger E and J, we expect the boson stars to eventually become unstable. For nonrotating 
boson stars this usually occurs at the maximum value of E, i.e. , at the first cusp. However, 
in the rotating case, it is clear from Fig. 4 that the maximum of E does not coincide with 
an extremum of the frequency u. Thus, if one looks for a zero mode to mark the transition 
between stability and instability, it will not occur at the maximum mass. 

We now turn to black holes. If we consider E and J in the region where MP-AdS black 
holes are unstable, then our hairy black holes are almost certainly stable within our m = 1 
scalar ansatz. They represent the endpoint of the superradiant instability for this mode and 
have larger entropy than the MP-AdS solutions. Outside this region of {E, J) the hairy black 
holes probably remain stable to m = 1 perturbations for a while, but then become unstable 
(like the boson stars). However, in all cases, since our hairy black holes have Qh^ > 4, it is 
clear from Fig. 1 that they are likely to be unstable to m > 1 perturbations. It is natural to 
ask what is the endpoint of this instability in the full theory. This raises an interesting puzzle. 

If the black hole settles down, the scalar fields at the horizon must satisfy i^^V^H = 
where is the null generator of the horizon. This ensures that there is no flux of scalar field 
into the black hole. If the final black hole has Qh^ > 1 it will probably still be unstable to 
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superradiance in a large m mode, so we need < 1- Since evolution is likely to stop at 
the first stable configuration, we expect the final black hole to have = 1. But since we 
can start with an arbitrarily small black hole, this would seem to require a boson star with 
Qh"^ = I5 i-G- invariant under dt + d^. Since this Killing vector does not have an ergoregion, 
such a boson star is unlikely to exist. 

The only alternative would be that the black hole will never settle down. It was shown in 



section (4.1.2) that the time scale for the (linear) instability grows rapidly with m. Thus, under 
evolution, the black hole approaches the m = 1 stable configuration discussed here. Then, the 
m = 2 instability kicks in and the black hole evolves toward the m = 2 stable configuration. 
This can repeat indefinitely, with the black hole developing structure on smaller and smaller 
scales over longer and longer times. However, since there is only a finite amount of energy that 
can be extracted from the black hole, one still expects the amplitude for the m*^ mode to go 
to zero, so it is not clear why it doesn't settle down. 

We have checked that standard MP-AdS black holes always have lower free energy than our 
hairy black holes. This is true for both the Helmholtz free energy F = E — TS (comparing 
black holes with the same temperature and angular momentum) and the Gibbs free energy 
G = E — TS — VLhJ (comparing black holes with the same temperature and angular velocity). 

Our ansatz for the scalar fields has an immediate generalization to any odd dimensional 
spacetime. We expect that the results discussed here will extend to higher dimensions. We 
also expect that most of the properties that we have found for black holes coupled to massless 
scalar fields will carry over to black holes with no scalar fields. Vacuum black holes exhibit 
superradiance and with AdS boundary conditions we expect them to evolve toward black holes 
with gravitational waves propagating around them. In this case, the metric itself (and not just 
the scalar hair) would be invariant under just a single Killing field. We cannot add a mass to 
our scalars and look for asymptotically flat flve dimensional black holes with scalar hair. There 
is a superradiant instability for massive scalars only in four dimensions [22]. 

We have considered only classical solutions, so it is natural to ask what happens if one 
includes Hawking radiation. Since our black holes are all small compared to the AdS radius, 
they have negative specific heat and one might expect them to evaporate completely. However 
this is incorrect since it would result in a decrease in entropy: A small spherical black hole with 
horizon radius r+ has energy ~ r^/G and entropy Sbh ~ r^/G. The entropy of a gas with 
energy E in AdS^ is roug hly Sg ~ {Eif^ where i is the AdS radius. If all the energy in the 
black hole went into the gas, the entropy would increase only if r+/L < (£p/£)^/^ where Ip is the 
Planck length. Our black holes can be much larger than this, so they will start to evaporate 
but then quickly come into equilibrium with their Hawking radiation. 

Given the recent work on the fiuid/gravity correspondence (see, e.g., PB]), one might wonder 
if there is a dual interpretation of our solutions in terms of a rotating fiuid. The answer appears 
to be no. The hydrodynamic approximation involves a restriction to long wavelength modes. 
Our black holes are all small compared to the AdS radius, so they do not describe a state of a 
fiuid. To describe a thermal state of a fiuid one needs a large black hole with r+ ^ i. Such a 
black hole can still have a superradiant instability, but only for modes with large m. This would 
again violate the long wavelength requirement of the hydrodynamic approximation. Even if 
there is no fiuid interpretation, there must still be some dual interpretation of our hairy black 
holes in terms of gauge/gravity duality, and it would be interesting to understand what it is. 



Note Added: After this paper was submitted, we were informed of [37] which describes 
another type of black hole with only one symmetry. It is static and asymptotically fiat, but 
has no spatial symmetries. 
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A Perturbative fields for the boson star and hairy black 
hole 

A.l Perturbative boson star 



In subsection 3.1 we explained how we can construct perturbatively the boson stars. In this 
Appendix we explicitly write the expansion of the gravitational and doublet scalar fields up to 
order O fe^): 



fir) 



9(^2 + ^2)4 1270080 (r2 + ^2)9 \^ 
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+0 {e') . (A.l) 
The perturbative construction also generates the expansion for the frequency of the scalar 



field, which is explicitly written in equation (3.4). 
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A. 2 Perturbative black hole: near-region and far-region field expan- 
sions 



In subsections |3.2.H|370] we have explained how we can construct perturbatively the hairy 
black holes using a matched asymptotic expansion of the near and far region solutions. In this 
Appendix we explicitly write the expansion of the gravitational and doublet scalar fields in the 
near-region and far-region with the matching conditions already imposed. We do the expansion 
up to order O (€"■, (r+/£)P) with {n + p) = 4. 

The far-region expansion of the gravitational and scalar fields is given by: 
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where Li„(a;) = ^^^iX^ /k"' is the Polylogarithm function. 

The near-region expansion of the gravitational and scalar fields is given by: 
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and E{x)= f\l-ify^'\l-xifY'^dr^, K{x) = - ,f) {l - xr]^)]'^^^ dr^ 

Jo Jo 

are, respectively, tlie complete elliptic integral and the complete elliptic integral of the first 

kind. 

The matched asymptotic expansion construction also generates the expansion for the fre- 
quency of the scalar field, that is equal to the angular velocity of the hairy black hole, and 
explicitly written in equation (3.14). 
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